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Oscillations of B° mesons have been studied in samples selected from about 3.5 
million hadronic Z decays detected by DELPHI between 1992 and 1995. One 
analysis uses events in the exclusive decay channels: B° — > D~7r + or D^a^ and 
Oh: B s ° -> WK- 7r + or D K &i, where the D decays are completely reconstructed. 

In addition, B° — B° oscillations have been studied in events with an exclusively 
reconstructed D s accompanied in the same hemisphere by a high momentum 
hadron of opposite charge. Combining the two analyses, a limit on the mass 
difference between the physical B° states has been obtained: 

Am B o > 4.0 ps" 1 at the 95% C.L. 
with a sensitivity of Am B o =3.2 ps -1 . 

Using the latter sample of events, the B° lifetime has been measured and an 
upper limit on the decay width difference between the two physical B° states 
has been obtained: 

r B o = 1.53j#i5(stat.) ± 0.07(syst.) ps 

Ar B o/r B o < 0.69 at the 95% C.L. 
The combination of these results with those obtained using D^£ T sample gives: 

Am B o > 4.9 ps" 1 at the 95% C.L. 
with a sensitivity of Am B o =8.7 ps -1 . 

r B o = 1.46 ± 0.11 ps and Ar B o/r B o < 0.45 at the 95% C.L. 



(Eur. Phys. J. C18(2000)229) 
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1 Introduction 

In this paper, the average lifetime of the B° meson is measured and limits are derived on 
the oscillation frequency of the Bg-B° system, Am B o, and on the decay width difference, 
AFbq, between mass eigenstates of this system. 

Starting with a B° meson produced at time t=0, the probability, V(t), to observe a B° 
or a B° decaying at the proper time t can be written, neglecting effects from CP violation: 

V(t) = ^|i e " rB e°'(cosh(^it) ± cos(Am B ot)) (1) 

where T B o = (r^o + T B , )/2, Ar B o = Tgo — Tg and Am B « = mj|o — mjj ; L and 
H denote the light and heavy physical states, respectively; Ar B o and Am B o are defined 
to be positive [Q and the plus (minus) signs in equation ([[]) refer to Bg (B°) decays. 
The oscillation period gives a direct measurement of the mass difference between the 
two physical states. The Standard Model predicts that APbo <C Am-QO and the previous 
expression simplifies to : 

V$T*(t) = rBoe-^^cos 2 !^^) (2) 
for Bg — > B° and similarly: 

\t) - L B oe B = sin ( — - — ) (6) 

for B^ — > B^. The oscillation frequency, proportional to Ameo, can be obtained from 
the fit of the time distributions given in relations (0) and (^), whereas expression (|Tj), 
without distinguishing between the Bg and the B°, can be used to determine the average 
lifetime and the difference between the lifetimes of the heavy and light mass eigenstatesQ. 



B physics allows a precise determination of some of the parameters of the CKM matrix. 
All the nine elements can be expressed in term of 4 parameters that are, in Wolfenstein 
parameterisation 0, A, A, p and rj. The most uncertain parameters are p and r\. 

Several quantities which depend on p and r\ can be measured and, if the Standard 
Model is correct, they must define compatible values for the two parameters inside mea- 
surement errors and theoretical uncertainties. These quantities are ex, the parameter 
introduced to measure CP violation in the K system, iKbl/lKbl, the ratio between the 
modulus of the CKM matrix elements corresponding to b^u and b— >c transitions, and 
the mass difference Am B o. 

In the Standard Model, B°-B° (q = d, s) mixing is a direct consequence of second 
order weak interactions. Having kept only the dominant top quark contribution, Attt-bo 
can be expressed in terms of Standard Model parameters : 

Am B o = ^\V th \ 2 \V tq \ 2 mlm B JlB BqVB F(x t ). (4) 

2 

In this expression Gp is the Fermi coupling constant; F(x t ), with x t = -^r, results from 

m w 

the evaluation of the second order weak "box" diagram responsible for the mixing and has 
a smooth dependence on x t ; rj B is a QCD correction factor obtained at next-to-leading 

1 Throughout the paper the dimension of Am B o will be expressed in ps -1 units, because the argument Am B ot in the 
expression (0) has no dimension. The conversion is given by ps _1 = 6.58 X 10~ 4 eV/c 2 . 
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order in perturbative QCD [||. The dominant uncertainties in equation (^) come from 
the evaluation of the B meson decay constant / Bq and of the "bag" parameter 5 Bq . The 
two elements of the Vckm matrix are equal to: 



\V td \ = A\^{l-pf + ^ , \V ts \ = AX 2 , (5) 

neglecting terms of order 0(X 4 ). 

In the Wolfenstein parameterisation, |V^ S | is independent of p and 77. A measurement 
of AmB s o is thus a way to measure the value of the non-perturbative QCD parameters. 

Direct information on V t d can be inferred by measuring Am B o. Several experiments 
have accurately measured Am B o , nevertheless this precision cannot be fully exploited to 
extract information on p and 77 because of the large uncertainty which originates in the 
evaluation of the non-perturbative QCD parameters. 

An efficient constraint is the ratio between the Standard Model expectations for Am B o 
and Am B o, given by: 

Am B o m B o/|o-B B o?7 B o \y |2 

°d _ d D d d d I r td| / ^? \ 

Am B o m B ofvoB B or] B o \V ts \ 2 ' 

A measurement of the ratio Am B o/Am B o gives the same type of constraint, in the p — rj 
plane, as a measurement of Am B o , but is expected to be better under control from theory, 

d 

since the ratios of the decay constants / Bq and of the bag parameters B Bq for and B° 
are better known than their individual values ||. 

Using existing measurements which constrain p and 77, except those on Am B o, the 
distribution for the expected values of Am B o can be obtained. It has been shown, for 
example, in the context of the Electroweak Standard Model and QCD assumptions, that 
Am B o has to lie, at 68% C.L., between 8 ps _1 and 16 ps -1 and is expected to be smaller 
than 21 ps" 1 at the 95% C.L. §]. 

The search for Bj? — B° oscillations has been the subject of recent intense activity. 
No signal has been observed so far and the lower limit on the oscillation frequency 
comes from the combination of the results obtained at LEP, CDF and SLD experiments: 
Am B o > 14.3 ps -1 at the 95% C.L. [0. The sensitivity of present measurements is es- 
timated as 14.7 ps -1 . These results have been obtained by combining analyses which 
select various B° decay channels^. Some analyses |||9| use events containing simply a 
lepton emitted at large transverse momentum relative to the axis of the jet from which 
it emerges. In this case, the proper time is measured using an inclusive vertex algorithm 
to reconstruct the decay distance and the energy of the B hadron candidate. In other 
analyses, like D^£ T [jT0| , pTTf1 , the identified lepton is accompanied, in the same hemisphere, 
by an exclusively reconstructed D s . In D^h^ analyses (see |12| and the present paper), 
such leptons are replaced by one or more charged hadrons. 

Progress before the next generation accelerators is expected to come from improved 
analyses of these channels, but the sensitivity at high frequency is essentially limited by 
the damping of the reconstructed oscillation amplitude due to the limited resolution on 
the B° proper time. In general, the proper time resolution a t contains two terms and one 
of them is a time- dependent: 



o t = s]°l 2 + (°p/p) 2 x t 2 (7) 

where <jl is related to the decay distance resolution, cr p /p is the relative error on the 
momentum reconstruction and t is the proper decay time. The damping factor of the 
oscillation amplitude in this case is given by the following expression 

2 Unless explicitly stated otherwise, charge conjugate states are always implied. 
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damping = exp(— (Am B o<r^ 2 /2)v / 7r)(x) exp(x 2 )ERFC(x), (8) 
where x=l/ v / 2/Am B o l/(<r p /p) and ERFC(a;)=2/ v ^F J^e^dt. 

Exclusive (completely reconstructed) decays have a better time resolution for two reasons. 
As there is no missing particle in the decay, the B|? momentum is known with a good 
precision and therefore the contribution of the momentum uncertainty to the proper 
time resolution is negligible. Hence, at = <Jl and the previous expression is simplified: 
damping = exp(— (Am B o<7L) 2 /2). In addition, the reconstructed channels are two-body 
or quasi two-body decays, with an opening angle of their decay products which is on 
average larger than in multi-body final states; this results in a better accuracy on their 
decay position determination. 

The B° meson lifetime is expected to be equal to the Bj lifetime [14| within one percent. 



In the Standard Model, the ratio between the decay width and the mass differences in the 
B°-B° system is of the order of (m^/mt) 2 , although large QCD corrections are expected. 
Explicit calculations to leading order in QCD correction, in the OPE formalism predict: 

Ar B o/r B o = o.i6i°;& 

where the quoted error is dominated by the uncertainty related to hadronic matrix ele- 



ments. Recent calculations [pl| at next-to-leading order predict a lower value: 



Ar B o/r B o = 0.006^° . 

An interesting approach consists in using the ratio between Ar B o and Am B o Hl5| : 



^ = (2.63^ 6 7 )KT 3 
Am B o 



(9) 



to constrain the upper part of the Am B o spectrum with an upper limit on Ar B o/r B o. If, 
in future, the theoretical uncertainty can be reduced, this method can give an alternative 
way of determining Am B o via Ar B o and, in conjunction with the determination of Am B o , 

s s d 

can provide an extra constraint on the p and n parameters. 

After a description of the main features of the DELPHI detector, the event selection 
and the event simulation in Section the two B° analyses are described. Section |^ 
presents an analysis of Bg-B° oscillations from a sample of 44 exclusively reconstructed 
B° decays. Section [| is dedicated to the D^h^ analysis, in which B^ — > D s hX decays 
with fully reconstructed D s are selected. This analysis also includes a measurement of 
the B° meson lifetime and sets a limit on the difference between the decay widths of 
the physical B^ states, Ar B o/r B o. In Section the combined limit on Am B o is given. 
Finally in Section || these three results are combined with those obtained using a 
sample flOfl. 



2 The DELPHI detector 

The events used in this analysis have been recorded with the DELPHI detector at LEP 
operating at energies close to the Z peak. The DELPHI detector and its performance 
were described in detail elsewhere ||16|| . In this section, components of the detector and 
their characteristics, which are the most relevant for these analyses are summarised. 
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2.1 Charged particle reconstruction 

The detector elements used for tracking are the Vertex Detector (VD), the Inner 
Detector (ID), the Time Projection Chamber (TPC) and the Outer Detector (OD). The 
VD provides the high precision needed near the primary vertex. 

For the data taken from 1991 to 1993, the VD consisted of three cylindrical layers of 
silicon detectors (at radii of 6.3, 9.0 and 10.9 cm) measuring points in the plane transverse 
to the beam direction (r0 coordinate^]) in the polar angle range between 43° and 137°. 
In 1994, the first and the third layers were equipped with detector modules with double 
sided readout, providing a single hit precision of 7.6 /im in r<fr, similar to that obtained 



previously, and 9 /im in z Wn - F° r high momentum particles with associated hits in the 
VD, the impact parameter precision close to the interaction region is 20 jum in the r<p 
plane and 34 /im in the rz plane. Charged particle tracks are reconstructed with 95% 
efficiency and with a momentum resolution a p /p < 1.5 x 10~ 3 p (p in GeV/ c) in the polar 
angle region between 25° and 155°. 

2.2 Hadronic Z selection and event topology 

Hadronic events from Z decays were selected by requiring a multiplicity of charged 
particles greater than four and a total energy of charged particles greater than 0.12y/s, 
where i/i is the centre-of-mass energy and all particles were assumed to be pions; charged 
particles were required to have a momentum greater than 0.4 GeV/ c and a polar angle 
between 20° and 160°. The overall trigger and selection efficiency is (95.0±0.1)% fi~8f . A 
total of 3.5 million hadronic events was obtained from the 1992-1995 data. 

Each selected event was divided into two hemispheres separated by the plane transverse 
to the sphericity axis. A clustering analysis, based on the JETSET algorithm LUCLUS 



with default parameters p9| , was used to define jets using both charged and neutral 
particles. These jets were used to compute the p° ut of each particle of the event, as the 
transverse momentum of this particle with respect to the axis of the jet it belonged to, 
after having removed this particle from its jet. 

2.3 Hadron identification 

Hadron identification relied on the RICH detector and on the specific ionisation mea- 
surement performed by the TPC. 

The RICH detector |16| used two radiators. A gas radiator separated kaons from 
pions between 3 and 9 GeV/c, where kaons gave no Cherenkov light whereas pions did, 
and between 9 and 16 GeV/c, using the measured Cherenkov angle. It also provided 
kaon/proton separation from 8 to 20 GeV/c. A liquid radiator, which has been fully 
operational since 1994, provided p/K/ir separation in the momentum range between 0.7 
and 8 GeV/c. 

The specific energy loss per unit length (dE/dx) was measured in the TPC by using 
up to 192 sense wires. At least 30 contributing measurements were required to compute 
the truncated mean. In the momentum range between 3 and 25 GeV/c, this is fulfilled 
for 55% of the tracks and the dE/dx measurement has a precision of ±7%. 

The combination of the two measurements, dE/dx and RICH, provides four levels 



of pion, kaon and proton tagging 20]: "very loose", "loose", "standard" and "tight" 



corresponding to different purities. The efficiency and purity of the tagging depend on 

3 The DELPHI reference frame is defined with z along the e — beam, x towards the centre of LEP and y upwards. 
Angular coordinates are 8, measured from z, and the azimuth, <j>, from the x-axis, while r is the distance from the z-axis. 
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the momentum of the particle. In the typical momentum range of the B decay products 
between 2 and 10 GeV/c, the average efficiency (purity) is about 75% (50%) for the 
"loose", 65% (60%), for the "standard" and 55% (70%) for the "tight" kaon tag. The 
"very loose" tag indicates that the particle is not identified as a pion. The purity is the 
proportion of genuine kaons in the selected sample. 

Very recently a neural network algorithm has been developed in order to combine the 
different RICH identification packages optimally. The result of the neural network gives 
a tagging variable x net , which varies from —1 (pure background) to +1 (pure signal). In 
order to use the same definitions of the kaon tag through all the paper, the "very loose" , 
"loose" , "standard" and "tight" tags were defined for a tagging variable x net larger than 
—0.2, 0.0, 0.4 and 0.6, respectively. This second algorithm has been used in the D^h T 
analysis (Section |). 



2.4 Primary and B vertices reconstruction 

The average beam spot position evaluated using the events from different periods of 
data taking has been used as a constraint in the determination of the e + e~ interaction 



point on an event-by-event basis [IE|. In 1994 and 1995 data, the position of the primary 
vertex transverse to the beam is determined with a precision of about 40 fim in the 
horizontal direction, and about 10 /im in the vertical direction. For 1992 and 1993 data, 
the uncertainties are larger by about 50%. 

For both 1992-1993 and 1994-1995 data, the B decay length was estimated as L = 
L xy / sin^B, where L xy is the measured distance between the primary vertex and the B 
decay vertex in the plane transverse to the beam direction and 9-q is the polar angle of 
the B flight direction, estimated from the B decay products. The transverse decay length 
L xy was given the same sign as the scalar product of the B momentum with the vector 
joining the primary to the secondary vertex; this procedure gives signed decay length L. 



2.5 A and K? reconstruction 



The A — > p7r~ and K° — > tt + tt~ decays have been reconstructed using a kinematic fit. 
The distance in the rcf) plane between the V° decay point and the primary vertex was 
required to be less than 90 cm. This condition meant that the decay products had track 
segments of at least 20 cm long in the TPC. The reconstruction of the V° vertex and 
selection cuts are described in detail in reference Hl6| . Only K° candidates passing the 



"tight" selection criteria defined in [16] were retained for this analysis. 



2.6 b-tagging 



The b-tagging package is described in reference The impact parameters of the 
charged particle tracks, with respect to the primary vertex, were used to build the prob- 
ability that all tracks come from this vertex. Due to the long B hadron lifetimes, the 
probability distribution is peaked near zero for events which contain beauty quark whereas 
it is flat for events containing only light quarks. 



2.7 Event simulation 



Simulated events were generated using the JETSET 7.3 program JT3 with parameters 
tuned as in P2| and with an updated model of B decay branching fractions. B hadron 
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semileptonic decays were simulated using the ISGW model p3 



Generated events were 
passed through the full simulation of the DELPHI detector [16], and the resulting sim- 
ulated raw data were processed through the same reconstruction and analysis programs 
as the real data. 



3 Exclusively reconstructed Bg decays 

3.1 Evaluation of Bg branching fractions 

For this analysis B° mesons were reconstructed in D~7r + , D^Ta^, D°K~7T + and D°K~a^ 
decay channels. Contributions to the mass spectrum can come also from other decay chan- 
nels such as D*~7r + , D*-a|, D*(2007)°K-7r + , D*(2007)°K-af , D"p + and D*"P + , where 
the cascade photon or the neutral pion has not been reconstructed. All correspond- 
ing branching fractions are unmeasured. An estimation will be used in the following to 
evaluate the number of the expected B° events. 

To evaluate the two-body B° — > D^ - n + (af) branching fractions, the equivalent decay 
channels for non-strange mesons were used. This evaluation is based on the numerical 



application of factorisation done in reference |24 . 

For each decay mode of the Bg — > DW~M + , where M + is a ix + or a p + , the branching 
fraction is estimated using the following formula: 

Br(B° -> DW-M+) = Br cxp (B° -> dW~M+) ' ) * * / , (10) 

v h ' y d Br th (B°j — > D(*)~M+) V ; 

where the notation D^ - means D*(2010)~ or D~ meson. The superscript "exp" des- 
ignates the experimentally measured values f25| , and "th" the predictions from refer- 
ence 



24f . The theoretical ratios in equation (|T0|) are close to 1. 
For the two-body decays containing D mesons with ai, the theoretical branching frac- 
tions are not available, and the same factors as for the p + channels have been assumed. 



Bjj decay channel 


Measured Br 


Bg decay channel 


Estimated Br 


Bo -> D-7T+ 


(3.0 ±0.4) x 1(T 3 


Bo _> D-7T+ 


(3.4 ±0.4) x 1(T 3 


B^ -> D*(2010)-7T+ 


(2.8 ±0.2) x 1(T 3 


Bo -> D*-7T+ 


(3.0 ±0.2) x 1(T 3 


Bo - D~p+ 


(7.9 ± 1.4) x 1(T 3 


Bo - D-p+ 


(9.0 ± 1.6) x 1(T 3 


B^ -> D*(2010)-p+ 


(6.7 ±3.3) x 1(T 3 


Bo - D s *-p+ 


(6.9 ±3.4) x 1(T 3 


Bo -> D-a^ 


(6.0 ±3.3) x 10~ 3 


Bo -> D s -a+ 


(6.8 ±3.6) x 10~ 3 


-> D*(2010)-a^ 


(13.0 ±2.7) x 1(T 3 


BO -> D s *-a+ 


(13.3 ±2.8) x 10~ 3 


B o ^ D W Vtt+ 


see Section |3~4| 


B o ^ d W0 K-7T+ 


(0.9 ±0.5) x 1(T 3 


B o _> D ( * ,0 vr-a+ 


see Section 3.4 


B« - D ( * )0 K-a+ 


(3.0 ±1.7) x 1(T 3 



Table 1: Values of the Bg branching fractions used in the current analysis. The notation 
defines D*(2007)° or D° charmed hadrons. 



Similar considerations are applied to the three-body decays Bg — > D K 7r + (a+) and 
B° -> D*(2007)°K-7r+(a^) using the measurements of Br(B° 1 -> D°7r-7r+) and Br(B° d -> 
D°7r _ a^), which are described in the Section |3.4| of the current paper. The ratio of Bg 
to B° branching fractions was taken equal to 1.1 (this factor corresponds to the average 
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ratio of the four previous decay modes). Finally, the factor 3.3 was used as the ratio 
between the B° decay modes containing af and ir + mesons. This factor was taken as an 
average between the experimental values for the following B° two-body decays: 



Br(B° 



D" 



Br(B° -> D*(2010)-a^ 



Br(B° 



D- 



7T 



Br(B° -> D*(2010)- 



3.3 ±0.9 



"7T 



The theoretical calculation gives for this ratio values between 3.4 and 3.5. Table [1| 
presents the evaluations of several branching fractions of interest for this paper. 



3.2 Event Sample 

Events were selected using the following decay channels: 

07T-, (f)TT-TT + n-, f(980)7r- K S °K~, K* K~, K*°K*- 

K + 7T~, K + 7T~7r + 7r~ 
K + 7T~, K + 7T~7r + 7T~ 

where the 0, K°, f(980), K*°, K*~ and ai are reconstructed in their charged decay channels: 
-> K+K-, K° -> 7T+7T-, f(980) -> vr+vr-, K*° -> K+tt", K*~ -> K°vr- and -> p°7r+, 
p° — > 7r + 7r _ . D s and D° mesons were reconstructed by considering charged particles in 
the same hemisphere with at least one VD hit. 

D meson candidates were accepted if their mass was within the intervals 1.93 — 
2.01 GeV/c 2 for D s and 1.83 - 1.90 GeV/c 2 for D°. The D decay length was required to 
be positive and the x 2 ~P r °bability of the fitted vertex to be larger than 10~ 5 (10~ 3 for 
the D° — > K7T7T7r decay mode). Different selection criteria were used for different decay 
channels according to the optimisations derived from dedicated simulated samples. They 
are described in the following: 



B s ° 
B s ° 
B s ° 
B? 



D-7T + 




D s &i 


D- 


D°K-7r+ 


D° 


D°K-a^ 


D° 



The (ft meson was reconstructed in the decay mode <fi — > K + K _ by taking all possible pairs 
of oppositely charged particles if at least one of them was identified as a "very loose" 
kaon. The invariant mass of these pairs had to be within ±12 MeV/c 2 of the nominal 
<p mass value The momenta of all three particles were required to be larger than 

1 GeV/c. In the decay of the D s meson into a vector (0) and a pseudoscalar meson (jr), 
helicity conservation requires that the angle ip, measured in the vector meson rest frame 
between the directions of its decay products and of the pseudoscalar meson, has a cos 2 ip 
distribution. The value of | cost/>| was required to be larger than 0.3. 

• D~ — > (/)TT~1T + 1T~ 



The (j) meson was reconstructed as in the previous channel. The momenta of all three 
pions had to be larger than 0.6 GeV/ c. 

•D s 7 -> f(980)vr- 



The f(980) meson was reconstructed in the decay mode f(980) — > 7i + 7i~ by taking all 
possible pairs of oppositely charged particles classified as pions. This channel suffers 
from a large combinatorial background which was reduced by selecting candidates having 
an invariant 7r + 7r~ mass within 15 MeV/c 2 of the nominal f(980) mass 12511 and a total 



momentum larger than 8 GeV/c. The momenta of all three pions had to be larger than 
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1 GeV/c. 



» D- -» K S °K- 

'0 



meson reconstruction has been described in Section |275|. In addition, the decay length 
of K° candidates had to be positive and their momentum to be larger than 2.5 GeV/ c. 
The momentum of the K~ candidate, identified as a "very loose" kaon, had to be larger 
than 1 GeV/c. 

• D- — > K*°K~ 

The K*° meson was reconstructed in the charged decay mode K*° — > K + 7r~. The K + 
candidate was required to be identified as a "very loose" kaon. The momenta of both 
particles had to be larger than 1 GeV/ c and the invariant mass of the pair had to be 



within ±40 MeV/ c of the nominal K mass [25|. The value of | cos^l (see the D s — > </>7r 



selection) had to be larger than 0.3. For the K~ candidate, the combinatorial background 
is higher than for the kaon coming from K*° resonance. The momentum of the K~ can- 
didate from D~ had to exceed 2.5 GeV/c and it had to be identified as a "loose" kaon. 

• D s 7 — > K*°K*~ 

The K*° meson was reconstructed as previously, but with an invariant mass within 
±60 MeV/c 2 of the nominal K*° mass |2f|. This mass interval was chosen larger than 
in the previous K*° selection, because of the additional constraint from the K*~ mass. 
The K*~ meson was reconstructed in the decay mode K* _ — > Kg7r~. The K° meson 
reconstruction was discussed previously. The invariant mass of K*~ candidates had to be 



within ±60 MeV/ c of the nominal K* mass |25 



> D° -> K+7T-, D° -f K + TT-TT + TT- 

The D° meson decay to K + 7r~ has been reconstructed by combining a kaon candidate, 
identified with the "loose" tag, with an oppositely charged pion with momentum larger 
than 1 GeV/c. The D meson decay to K + 7r~7r + 7r _ has been reconstructed by com- 
bining a kaon candidate, identified with the "standard" tag, with three pions, each of 
them having a momentum larger than 0.5 GeV/c. In order to reduce the combinatorial 
background, for both decay modes, the kaon momentum was required to be larger than 
2.5 GeV/ c and the D momentum candidate to be larger than 10 GeV/ c. 

Selected D° and D~ mesons were used to reconstruct Bg candidates by fitting a com- 
mon vertex for the D~7r + , D~ai , D K~7r + or D K~ai~ systems. The 7r momentum had 
to be larger than 4 GeV/c. The ai candidates were reconstructed using the combination 
of three pions with momenta larger than 0.8 GeV/ c and with an invariant mass situated 
within the interval 0.95-1.50 GeV/c 2 . At least one of the two 7r + 7r~ combinations was 



required to have an invariant mass lying within ±150 MeV/c of the nominal p mass |25]. 
The ai momentum had to be larger than 5 GeV/c (6 GeV/c for the D° — > K + 7r~7r + 7r + 
channel). For all candidates, the D° and D~ meson decay distance, relative to the B 
vertex had to be positive. Events with an estimated error on the Bg decay distance 
larger than 250 \im and those having a vertex x 2 -probability smaller than 10 -3 were 
removed. In order to reduce the combinatorial background from charm and light quarks, 
the b-tagging probability for the whole event and for the hemisphere opposite to the 
reconstructed B meson had to be smaller than 0.1. Additional selections which depend 
on the Bg decay channel were applied, mainly for D decays which suffer from a larger 
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combinatorial background than D s candidates: 



• B»->D-7r+, B s "-.D-a+ 
The momentum of the B° had to be larger than 22 GeV/ c. For B° — > D~a^ candidates, 
only the combination with the largest B° momentum has been kept. 

•B? -> D°K-tt + , B° -> D°K-at 



For B° — > D°K~7r + decays, the B° momentum had to be larger than 27 GeV/ c. Because 
of a high combinatorial background for B° — > D°K~ai~ decays, the B° momentum had to 
be larger than 29 GeV/c in the D — > K + 7r~ channel and larger than 33 GeV/c in the 
D° — > K + 7r _ 7r + 7r _ decay channel. In each event, only the candidate with the largest B° 
momentum was kept. For the K~ candidate, identified as "standard", the momentum 
had to be larger than 2 GeV/c. The source of the D and K~ meson pair can be an 
excited D meson state. The constraint on the mass value of such a state allows more 
soft K~ mesons to be considered. The selection on the K~ momentum was reduced 
to 1.5 GeV/c, if the invariant mass of the D K~ pair was compatible with the mass 
of the orbitally excited D s j(2573) state (the most probable value of J is 2 |25|| ). One 
event was detected, which is compatible with the decay channel B° — > D s j(2573)~7r + , 

D s j(2573)~ — ► D°K~: the momentum of the K~ candidate is 1.8 GeV/c and the mass 
difference M(D°K") - M(D°) is 716.0 ± 2.1 MeV/c 2 . The expected mass difference of 
708.9 ± 1.8 MeV/ c 2 is in good agreement with the observed one, taking into account the 
full width of this state which is 15±| MeV/c 2 J25 . 



Main peak 


Data set 


Mass 
(GeV/c 2 ) 


Width (a) 

(GeV/c 2 ) 


^signal 


Comb, bkg 

% 


Reflection 

% 


Real Data 


5.373 ±0.016 


0.029 ±0.012 


8±4 


27 ± 16 




Simulation 


5.370 ±0.002 


0.037 ±0.002 


5± 1 




7±3 


Satellite peak 


Real Data 


5.050 ±0.054 


0.111 ±0.049 


15 ±8 


55 ± 13 




Simulation 


5.099 ±0.007 


0.148 ±0.007 


11 ±2 




12 ±6 



Table 2: Characteristics of the B° signals and comparison with the simulation. The 
expected numbers of events were calculated taking into account the different branching 
fractions as given in Table [1 and the corresponding reconstruction efficiencies. These 
efficiencies vary from (0.2 ± 0.1)% for the B° — ► D K^a^ with D — > K + n^n + n^ , up to 
(10.2 ± 0.8)% for the B° — > 0^^+ with D^ — > (p-n~ . The fitted or expected number of 
signal events and the fraction of combinatorial background are given inside a mass window 
corresponding to ±3o~(±2o~) for the main (satellite) peaks. In real data the number of 
events in these mass windows are 11(33). The B° mass in the simulation is 5.370 GeV/ c 2 . 
In the last column the fraction of physical background, discussed in detail in Section \3.!% 
is given. 



The mass spectrum, obtained by summing up the contributions from the different 
channels is shown in Figure |l|. The data are indicated by points with error bars and the 
fit result is shown by the solid line. The mass distribution in the signal region was fitted 
using two Gaussian functions of different widths to account for the exclusive B° signal 
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(main narrow peak) and for the presence of partially reconstructed B° decays (satellite 
wider peak). According to the expected branching fractions, it was assumed that the 
dominant contributions for the satellite peak come from the following decay channels: 
D*-7r+, D*-a|, D*(2007)°K-7r+ D*(2007)°K- a ^, D s -p+ and D*"p + , where the cascade 
photon or neutral pion (or both) were not reconstructed. In the simulation, it was verified 
that the mass distribution of the satellite peak can be described by a Gaussian function 
as shown in Figure |2|. The central values and the widths of the Gaussian functions 
(main and satellite peaks) in Figure [I] were left free to vary in the fit. The combinatorial 
background was fitted using an exponential function with a slope fixed according to the 
simulation. This slope was verified in several ways. Using 2.9 million bb and 5.3 million 
qq simulated events, the expected mass spectrum was obtained after having removed the 
B° signal contribution in the mass region of the satellite and main peaks. The sum of the 
contributions from b events, charm events and light-quark events is shown in Figure [l| 
and is in agreement with that obtained in data. Two further checks were performed, 
using events selected in the side-bands of the D~ and D candidates and from wrong sign 
combinations. The slopes of both distributions are in agreement with those obtained by 
fitting the data and the simulation. 

The fit to the mass distribution yielded a signal of 8 ±4 B° decays in the main peak and 
15 ± 8 events in the satellite peak. The probability that the background has fluctuated 
to give the observed number of events in the main peak is 3 x 10~ 4 . Table ||] gives 
the characteristics of the observed signals and the comparison with simulation. Specific 
decay channels from B°, B°, B + and A° were simulated with statistics ranging from ten 
to several thousand times the expected rates in real data. These samples were used to 
determine efficiencies and reflection probabilities discussed below. 

3.3 Reflections from B[j and A£ decays 

For several B° decay channels a possible physical background originates from non- 
strange B decays, when one of the pions or proton is misidentified as a kaon (kinematic 
reflections). The main decay channels are B[J — > D _ 7r + (a ] f ), D~ — > K*°7r~ (which can 
contribute to B° — > D ! 77r + (a ] f ), D~ — > K*°K~ candidates) and from B^ — > D°7r _ 7r + , 
D° — > K + 7r~ or K + 7r~7r + 7r~ (which can contribute to B° — > D°K - 7r + candidates). In 
the mass region of the main peak, 0.32 ± 0.13 events are expected to originate from 
kinematic reflections (with a 19% contribution from A° decays). These estimates were 
obtained using dedicated Monte Carlo samples of B[| and A° events passed through the 
full reconstruction algorithms and satisfying the selection criteria. The corresponding 
uncertainties come from the limited knowledge of the assumed branching fractions and 
from the simulated events statistics. 

To study the contribution in the satellite peak from kinematic reflections, the same 
decay processes were considered with channels in which the D meson is accompanied by 
a p. This gives 1.3 ± 0.7 events. 

In order to look for possible signals coming from kinematic reflections in real data, 
B° candidates were considered in turn as B° or A° hadrons by changing the kaon into a 
pion or an antiproton, respectively. The mass distributions obtained are similar to those 
expected from genuine simulated B° mesons and do not show any accumulation of events 
in the B° or A° mass regions. 
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3.4 Reconstruction of non-strange B meson decays 

Non-strange B mesons, decaying into a D° and a small number of pions, were re- 
constructed in order to verify the B° reconstruction algorithms. The B + — ► D 7r + , 
B d — ► D*(2010) _ 7r + and B d — > D*(2010)~a ] f decay channels were studied and their mass 
distributions are shown in Figures |3|a and [3]b. They were fitted using a Gaussian func- 
tion for the signal and using an exponential function for the combinatorial background. 
An additional Gaussian function was used in Figure 0a to account for the signal coming 
from the following decay channels^ B+ -> D*(2007)°7r+, B+ -> D°p+, B+ -> D*(2007)°p + , 
where the ir° and/or 7 from the D (2007)° and/or p decays were not reconstructed. The 
main selection criteria which were imposed are rather similar to those used in the B° 
analysis. The information relevant to these reconstructed channels is given in Table |3|. 
The numbers of observed events are in agreement with expectations. 



Channel 


B meson mass 
(GeV/c 2 ) 


N obs 


N 


B + _> dV 


5.278 ±0.008 


24 ±5 


28 ±5 


B d — > D*(2010)~7r + + D*(2010)~a^ 


5.277 ±0.011 


11 ±5 


10 ±4 


B° -> D U 7T-7r+ + DVat 


5.287 ±0.024 


8±4 





Table 3: Characteristics of the non-strange B decay mode and comparison with the 
simulation. The error on the expected number of events comes from the errors in the 
branching fractions and reconstruction efficiencies. The branching fraction of the decay 
mode in the last row is estimated in this paper. 

Finally two B° decay channels, B° A — > D°7r~7r + and B d — > D°7r _ a^ which are not yet 
well established, were considered. 

In the strange B sector (Section |3.1|) they correspond to the B° — > D K~ir + and 

B° -> D K~af decays. mesons were reconstructed using the same selection criteria as 
for the corresponding B° decay channel, but replacing a K meson by a 7r. All particles, 
not explicitly identified as protons or kaons were accepted as pions. In order to remove 
the B° d -> D*(2010)-vr+ and B° d -> D*(2010)"a^ contamination, candidates were required 
to have a mass difference M(D°7r~) — M(D°) larger than 0.16 GeV/c 2 . 

Figure [| shows the mass distribution for these two decay channels. It was fitted with 
an exponential function for the background and two Gaussian functions for the signals 
corresponding to the main and satellite peaks, with respective widths equal to about 
34 MeV/ c 2 and 60 MeV/ c 2 , as obtained from the simulation. The position of the wider 
Gaussian, corresponding to B° — > D*(20 07) 7r~7r + and B° d -> D (2007)° 7r a^ signals with 
unreconstructed ir° or 7, and the slope of the background exponential function were also 
taken from the simulation. The fit to the mass distribution yielded a signal of 8 ± 4 B d 
decays in the main peak. The probability that the background has fluctuated to give the 
observed number of events is 3 x 10~ 2 . 

The number of observed events in the main peak can be translated into a branching 
fraction: Br(B[J — > D 7t~tt + + D 7r~a^) = (3.6 ± 1.9) x 10~ 3 . The quoted error is com- 
pletely dominated by statistics. The CLEO Collaboration has given an upper limit p5| 
only on the decay channel with a pion pair: Br(B° — > D°7r~7r + ) < 1.6 x 10~ 3 (at the 
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7o C.L.). This limit is consistent with those which can be extracted from the measured 
channels in this paper: 

Br(B° d -> D°tt-tt + ) < 1.6 x 1(T 3 at the 90% C.L. 
Br(B° -> D 7r- ai + ) < 5.4 x 10~ 3 at the 90% C.L. 

For this evaluation it was assumed (see Section |3.1| ) that the decay channel with a^ is 
produced with a rate which is (3.3 ± 0.9) times larger than the channel with 7r + . 

3.5 Study of B°-Bj oscillations 

3.5.1 Algorithm for tagging b or b quark at production time 

The signature of the initial production of a b (b) quark in the jet containing the B° 
or B° candidate was determined using a combination of different variables sensitive to 
the initial quark state. For each individual variable X iy the probability density func- 
tions /bpQ) (/b(-^i)) f° r b (b) quarks were built and the ratio Ri = f^(Xi) / f h (Xi) was 
computed. The combined tagging variable is defined as: 

^ jj> 

x tag = -i : p , where R = T[Ri. (11) 
1 + ri 

The variable x tag varies between -1 and 1. Large and positive values of x tag correspond to 
a high probability that the produced quark is a b rather than a b in a given hemisphere. 

A set of nine discriminant variables has been selected for this analysis. One set (three 
variables) is determined in the hemisphere which contains the B° meson, the other set 
(five variables) in the hemisphere opposite to the B° meson, and one variable is common 
to both hemispheres. Details concerning the definition of these variables and the method 



are given in reference [10 



An event was classified as a mixed or as an unmixed candidate according to the sign, 
Qd-i of the D s electric charge for decay channels containing a D s meson, or according 
to the sign of the kaon for D K~7r + decay channels, relative to the sign of the x tag 
variable. Mixed candidates are defined by requiring Xtag-Qo < 0, and unmixed candidates 
correspond to x tag • Qd > 0. The probability, e^ g , for tagging correctly the b or b quark, 
from the measurement of x tag , was evaluated by using a dedicated simulated event sample 
and was found to be (74.5±0.5)%. 

The corresponding probability for events in the combinatorial background was ob- 
tained using real data candidates selected in the side-bands of the D signal: the probabil- 
ities to classify these events as mixed or unmixed candidates are called e™^ 6 and e"^ b iX , 
respectively. For the reflection events the analogous probabilities are called e™?, e!^p ix 
and their values were taken from the simulation. 



3.5.2 Proper time resolution 

For each event, the B° proper decay time was obtained from the measured decay length 
(Lbo) and the estimate of the B° momentum (pb°)- The measured position of the D~ 

or D° decay vertex, the momentum, and the corresponding measurement errors, were 
used to reconstruct a D~ or a D particle. A candidate B° decay vertex was obtained by 
intercepting the trajectory of this particle with the other charged particle tracks which 
are supposed to come from the B° decay vertex. 

For the main peak, the B° momentum is precisely known since all decay products are 
reconstructed. For the satellite peak, this is no longer true, because there are one or two 
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undetected neutral (s) (n° and/or 7) in the B° decay. As discussed in Section |3T2| , the satel- 
lite peak is assumed to be composed of D* _ 7r + , D*~a+, D~p + , D* _ p + and D*(2007)°K~7r + , 
D (2007)°K _ af decays. The first four decay channels give D~7r + (a^) and the last two 

channels give D°K _ 7r + (a ] f ) final states. Kinematic mass constrained fits (imposing the 
mass of B° and of the intermediate states with corresponding errors) were performed 
assuming always the presence of a 7 coming from the D*~ or D (2007)° meson decays. 
This is a good approximation also for decay channels containing a p + meson. These mass 
constrained fits were performed on events lying in the mass region corresponding to ±2<r 
of the fitted mass of the satellite peak. 

Except for the combinatorial background contribution, the predicted proper time dis- 
tributions were obtained by convoluting the theoretical functions with resolution functions 
evaluated from simulated events. Due to the different decay length resolutions (for the 
different Vertex Detector configurations), proper time resolutions were used for data sam- 
ples taken in 1992-1993 and 1994-1995 separately. The proper time resolution, TZB s {t — ti), 
was evaluated from the difference between the generated time (t) and the reconstructed 
time (tj), fitting this distribution by the sum of two Gaussian functions: 

n Bs (t-U) = {I- f 2 )G l {t-t l ,a 1 ) + f 2 G 2 {t-t h a 2 ) (12) 

where G\(t — ti,ai) and G 2 (t — ti,a 2 ) are Gaussian functions with the corresponding 
resolutions o\ and a 2 . In the general case, the Oi resolutions depend on the momentum 
uncertainty (see equation (^)). But in this analysis, as the B° momentum is known 
with a good precision, the contribution of the momentum uncertainty to the proper 
time resolution is negligible. f 2 is the fraction of the second Gaussian function, which, by 
convention, is that with the worse resolution. The values of the corresponding parameters, 
obtained from simulated events, are given in Table |j. The time resolution for events in 
the satellite peak is only slightly worse than for those belonging to the main peak. 



Main peak 


Decay channels 


Years 


cxi(ps) 


a 2 (ps) 


h 


all Bg channels 


1992-1993 


0.068 


0.18 


0.27 


all Bg channels 


1994-1995 


0.065 


0.12 


0.42 


Satellite peak 


D s *-7r+(af), D >P (2007) K-7r+(a^) 


1992-1993 


0.066 


0.15 


0.48 


D s *-7r+(a+), D >P (2007) K-7r+(a+) 


1994-1995 


0.081 


0.17 


0.30 


D-p+, B*-p + 


1992-1993 


0.085 


0.21 


0.65 


d- p + , d:- p + 


1994-1995 


0.092 


0.20 


0.57 



Table 4: Proper time resolution from simulation for 1992-1993 and 1994-1995 data 
sets and for the different decay channels. It has been parameterised by the sum of two 
Gaussian functions. 

The time distribution V CO mb{ti) for the combinatorial background under the main peak 
is obtained from real data by fitting the time distribution of wrong-sign and right-sign 
events lying in the side-bands of the B^ mass distribution. It was verified, in the simula- 
tion, that the time distribution for these classes of events is similar to that obtained for 
events lying under the B° mass peak. The time distribution V CO mb{ti) for the combina- 
torial background under the satellite peak was taken directly from the simulation, since 
it has a dependence on the measured Bg mass, due to the procedure used to reconstruct 
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the B momentum for these events. The time distribution V re f{ti) for the reflection was 
also taken from the simulation. 



3.5.3 Fitting procedure 



The oscillation analysis is performed in the framework of the amplitude method [13 



which consists in measuring, for each value of the frequency Am B o, an amplitude A and 



T=rO 



s 



its error cr{A). The parameter A is introduced in the time evolution of pure B° or B 
states, so that the value A = 1 corresponds to a genuine signal for oscillation. The 
time-dependent probability that B° is detected as a B° or B° is then: 

V unmix(mix)^ = _L x exp(-t/r B o) x (1 ± A cos(Am B ot)) (13) 
2r B o 

where the plus (minus) signs refer to B° (B°) decays. The 95% C.L. excluded region for 
Am B o is obtained by evaluating the probability that, in at most 5% of the real 
signal having an amplitude equal to unity would give an observed amplitude smaller than 
the one measured. This corresponds to the condition: 

A(Am B o) + 1.645 a(A(Am B o)) < 1. 

In the amplitude approach, the error cx(A(Am B o)) is related to the probability to exclude 
a given value of Am B o. The sensitivity is defined as the value of Am B o which would just 
be excluded, if the value of A, as measured in the experiment, was zero for all values 
of Am B o, i.e. it is the expected 95% confidence limit that an experiment of the same 
statistics and resolution would be expected to achieve (a real experiment would set a 
higher or lower limit, because of fluctuations in the values of A as a function of Am B o). 
The probability distributions for mixed and unmixed?] events are: 

P mX (U) = faPg'iU) + frefP™T(ti) + f COmb P™ b {U) (14) 

where / Bs , f re f and f CO mb are the relative fractions of the B°, reflection and combinatorial 
background events, respectively, which satisfy the condition / Bs +/ re /+/ C omb = l- The 
expressions for the different probability densities are: 

• B° mixing probability: 

Pb* x (U) ={ e b a9 V^ x (t) + (l-el a9 )V^ mx (t) }®ft Bs (t-t,,) (15) 

• Reflection mixing probability: 

PmU) = e?:?V ref (U) (16) 

• Combinatorial background mixing probability: 

Pcombi^i) = e combPcomb(ti) (17) 



The parameters e^ 9 , e™y and e™^ 6 were defined in Section |3.5.1 



The amplitude analysis is then performed using all events selected in a mass region 
between 4.83 and 5.46 GeV/ c 2 . The variation of the background level as a function of 
the reconstructed mass was included in this analysis on an event-by-event basis. Figure [5| 
shows the variation of the measured amplitude as a function of Am B o. With the present 

4 In the following, only the probability distribution for mixed events is written explicitly; the corresponding probability 
for unmixed events can be obtained by changing e into (1 — e). 
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level of the statistics, this analysis provides a negligible low limit. On the other hand, 
its most important feature, despite the low statistics, is the relatively small error on the 
amplitude at high values of Awib" with respect to the more inclusive analyses. Due to 
the limited statistics, the error on A (a (A)) can be asymmetric. The error has been 
symmetrized by taking the larger value. Figure |6|a shows the variation of the measured 
error on the amplitude as a function of Arneo. In Figure ^|b the results from the T)f£ T JTD 



and from the D^h T (see the second part of the current paper) are compared with the 
exclusive B° analysis. It should be noted that, due to the better proper time resolution, 
the resolution u{A) of the B° exclusive analysis increases more slowly. The ratio of the 
corresponding errors of B° exclusive and ~Df£ T analyses is about 5 (2) at the low (high) 
values of Ameo. 

The behaviour of <j(A) was investigated using a "toy" Monte Carlo generated with 
the same characteristics as those measured in real data. The individual toy experiments 
show a similar behaviour of cr(A), as in the real data, as a function of Am-go. For each 
value of Am-BO, the distribution of cr(A) and its variance was obtained. The central value 
and the region corresponding to a ±2a(A) variation are also shown in Figure |^a. 

3.5.4 Study of systematic uncertainties 

Systematic uncertainties were evaluated by varying the parameters, which were kept 
constant in the fit according to their measured or expected errors. 

• Systematics from the tagging purity: 
a variation of ±3% on the expected tagging purity for the signal was considered, 
following the results given in |10| . 



Systematics from the background level and kinematic reflection: 

the levels of background and of kinematic reflection were varied separately for the 

main peak and for the satellite peak according to the statistical uncertainties given 

in Table |2|. Since the level of the combinatorial background was used as a function 

of the reconstructed mass on an event-by-event basis, the measured central mass 

positions and the corresponding widths were also varied by ±1<t around their fitted 

values. 

Systematics from the expected resolution on the B decay proper time: 

the widths <j\ and <T2 described in Section 3.5.2 and given in Table [| were varied by 



±10% Hi3Ii 



• Finally, uncertainties in the values of the B° meson lifetime (tbo = 1.46±0.06 ps 0) 
and in the parameterisation of the proper time of the combinatorial background were 
taken into account. 

In Figure |5], which presents the variation of the measured amplitude as a function of 
AmBO, the shaded area shows the contribution from systematics. 

4 D^h T analysis with fully reconstructed D s 

Events with an exclusively reconstructed D s accompanied by one or more large momen- 
tum hadrons were used to perform a second oscillation analysis. This channel is similar 
to the final state ]TU[ but, instead of a charged lepton, it uses charged hadrons. 



It provides a larger number of events but suffers from an ambiguity in the choice of the 
hadrons and from a reduced B° purity of the selected sample. This approach has already 
been used in DELPHI to measure the B? lifetime lETi 



1G 



4.1 Event selection 

The D s meson is selected in two decay channels: 

D b 7^07r-, 0^K+IC 



Dr -> K* K~ K*° -> K 



s 



7T 



D s candidates were reconstructed by making combinations of three charged particles 
located in the same event hemisphere, with a momentum larger than 1 GeV/ c and with 
reconstructed tracks associated to at least one VD hit. The value of |cos^| (see Sec- 
tion |3.2|) was required to be larger than 0.4 and 0.6 in the D~ — > 07r~ and D~ — > K*°K~ 
decay channels, respectively. In order to reduce the combinatorial background from charm 
and light quarks, the b-tagging probability for the whole event was required to be smaller 
than 0.5. 

In this analysis the neural network algorithm for hadron identification described in 
Section |2]^ was used. For the D~ — > <p7r~ decay mode, the invariant mass of <ft candidates 
was required to be within ±12 MeV/ c 2 of the nominal <fi mass and the <fi momentum was 
required to be larger than 5 GeV/ c. If the K + K~ invariant mass was within ±4 MeV/ c 2 
of the nominal (f> mass, both kaon candidates were required to be identified as "very loose" 
kaons, otherwise to be identified as "loose" kaons. 

For the D~ — > K*°K~ decay mode, the invariant mass of the K*° candidates was re- 
quired to be within ±60 MeV/ c 2 of the nominal K*° mass value and the K*° momentum 
was required to exceed 5 GeV/ c. The momentum of the K~ candidate from D~ was re- 
quired to exceed 3 GeV/ c. To suppress the physical background from the D~ — > K + 7r~7r~ 
kinematic reflection, the K~ candidate was required to be identified as a "standard" kaon. 
For K*° — > K + 7r~ decays, "loose" identified K + were also accepted. In order to suppress 
the combinatorial background, the K~ candidate was required to be identified as "tight" 
kaon, if the invariant mass of the K*° candidates was out of ±20 MeV/ c 2 of the nominal 
K*° mass value, and the value of | cosip\ was required to be larger than 0.8. 

The selection of a hadron accompanying the D s candidate is based on an impact 
parameter technique. A sample of tracks coming predominantly from B hadron decays 
was preselected by using their impact parameters and the corresponding errors, both 
with respect to the primary vertex (Im p , a p ) and to the secondary D s decay vertex 
(Im s , a s ). The hadron was then searched for amongst the preselected particles in the 
event, by requiring that its charge was opposite to the D s charge and that it had the 
largest momentum. The efficiency of the hadron selection was about 80% and, among 
the selected hadrons, (84 ±4)% came from a B vertex. Details on the track preselection 
as well as on the hadron selection are given in reference |27 . 



The B decay vertex was reconstructed by fitting the selected hadron and the D s can- 
didate to a common vertex. The x 2 ~P r obability of the fitted B° vertex has been required 
to be larger than 10~ 3 . In order to increase the resolution on the measured decay length, 
only reconstructed events with a decay length error smaller than 0.07 cm were kept. 

If the previous procedure failed, a new attempt was made, using an inclusive algorithm 
which allowed several hadrons to be attached to the D s candidate. This algorithm is based 
on the difference in the rapidity distributions for particles coming from fragmentation 
and from B decays. The fragmentation particles, on average, have lower rapidity 



than B decay products. The charged particles were ordered in increasing values of the 
rapidity and were attached in turn to the secondary D s vertex. Up to three particles 
were accepted with their total charge equal to ±1 or 0. The rapidity was calculated 
as 0.5 log ((E + Pl)/(E — Pl)), where E is the energy of the particle (assumed to be a 
pion) and Pl its longitudinal momentum with respect to the thrust axis of the event. 
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Only particles with a momentum greater than 1 GeV/c were accepted. In addition, for 
the tracks satisfying the condition Im p /a p < 3, it was required that Im s /a s < Im p /a p . 
Events with a decay length error smaller than 0.07 cm and a ^-probability of the fitted 
B° vertex larger than 10~ 3 were kept. 



D s decay channels 


D s signal in 1992-1993 data 


D s signal in 1994-1995 data 




322 ±30 (0.60 ±0.04) 


468 ± 42 (0.53 ± 0.04) 


D" -> K* U K" 


152 ± 28 (0.70 ±0.06) 


324 ±35 (0.58 ±0.05) 



Table 5: Number o/D s mesons reconstructed in the 07r~ and K*°K~ decay channels after 
selection of the accompanying hadron(s). The fraction of combinatorial background, given 
in parentheses, has been evaluated using a mass interval of ±2<r centred on the fitted D s 
mass. 

The selected sample, specified as D^h^ in the following, contained about 30% of such 
"multi-hadron" events. Among them, about 20% have one, 50% two and 30% three se- 
lected hadrons. The multi-hadron and single-hadron events were treated in a similar way. 
Figure [7] shows the D s signals after selection of the accompanying hadron(s). The mass 
distributions were fitted with two Gaussian functions of equal widths to account for the 
D~ and D signals and with an exponential function for the combinatorial background. 
All parameters were allowed to vary in the fit. Table |5] gives the number of observed 
events in the D s signal, after background subtraction, and the fraction of combinatorial 
background. The D s signal region corresponds to a mass interval of ±2<r centred on the 
fitted D s mass. 

The reconstructed number of D~ — > 07r~ events using 1994-1995 data is about 1.5 
larger than those obtained using 1992-1993 data. This factor reflects the difference in 
statistics between the two data sets. Such "statistical scaling" does not apply to the 
D~ — > K*°K~ decay mode, where the particle identification, which was better for the 
1994-1995 data set (see Section |2|), plays a more important role. 

Four events are in common with the exclusively reconstructed B° sample: one event 
in the main peak and three in the satellite peak. These events were removed from the 
D^h^ sample for the oscillation analysis. 

4.2 Sample composition 

The D^h T sample contained a large physical background due to D s from non-strange 
B hadron decays and from cc fragmentation. Four sources of events, originating from B 
decays, were considered: two from B° and two from non-B^ mesons, namely, B decaying 
to one or two charmed mesons comprising at least one D s . The relative fractions of these 
sources were calculated using five input parameters: 

• Br(b -> DfX) at LEP 

• Br(b -> DjX) at T(4S) [Q; 

• Br(b -> B°) at LEP ||; 

• the probability in the non-strange B meson, that a D~ is produced at the lower 
vertex: Br (B u , d ^D-X) @; 

• the probability to have two charmed hadrons in a b-decay: Br(b ^DDX) 

The last two probabilities were assumed to be the same for all B species. To estimate 
the first two branching fractions , the averaged production rate of D s from all B species, 
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Br(b -> DfX) xBr(D s ± -> 0tt ± ), measured by the ALEPH, DELPHI and OPAL collab- 
orations and the equivalent quantity Br(B ud — > D^X) x Br(D^ — >■ (j)^), measured 
at the T(4S) by the CLEO and ARGUS collaborations pi| , were used. The following 
fractions for the different B decays contributing to the D^h T signal were evaluated: 

• Fraction of B° decaying to a D s and no other charmed meson: -Fb s ,id=(39 ± 7)%. 

• Fraction of B° decaying to a D s and another charmed meson: -Fb s ,2D=(H ± 3)%. 

• Fraction of B mesons (non-B^) decaying to a D s and no other charmed meson: 
F b ,id=(9±5)%. 

• Fraction of B mesons (non-B^) decaying to a D s and another charmed meson: 
F B ,2D=(41±7)%. 

The contribution F cc from direct charm was estimated from the measurement of D q 



production in charm events at LEP |30|, taking into account the Z partial widths into b 
and c quark pairs: F cc =(27± 5)%. 

Finally, the relative proportion of the combinatorial background f^g was taken directly 
from the fit of the real data mass distributions (see Table §). 



4.3 Discriminant variables to increase the Bg purity 

To increase the effective purity in B° of the selected sample, five variables were used 
which allow the separation of the signal and background components. These variables 
are: the D s mass, the D s momentum, the value of | cos^|, the ^-probability of the fitted 
D s decay vertex and the value of the b-tagging variable measured in the hemisphere 
opposite to that of the D s meson. 

The relative components in the selected sample, defined in the previous section, were 
calculated on an event-by-event basis: 

^b b ,id(1 - fee - /bfcg^bO 3 - tag) n T hc {Vi)/t0t 

i=l,4 

^b s ,2d(1 - fcc - fbkg)Fb$ ~ tag) n F hc {Vi)/t0t 

i=l,4 

-Fb,id(1 - fcc - fbk g )^b(b - tag) Y[ Fbc(vi)/tot 

i=l,4 

^b,2d(1 - fcc - fbkg)Fh{^ - tag) Y[ Fh C {Vi)/t0t 

i=l,4 

/cc^c(b-tag) J| F hc (vi)/tot 

i=l,4 

fbkg II Fbkg{Vi)/t0t 
i=l,5 

where Vi indicates the z-th discriminant variable, jF bc , jF b , T c and Tbkg are the probability 
density functions for the b and c together, b, c and the combinatorial background events, 
respectively. The relative charm contribution is f cc = F cc (l — fbkg)- In these expressions, 
the total normalisation factor is: 

, , feff I feff I feff I feff I feff I feff 

101 — JB Si m ' Jb s ,2D ~t Jb,1D + Jb,2D + ice + Jbkg ■ 

All discriminant variables, except for the b-tagging, were used to separate b and c events 
together from the combinatorial background (bkg). The b-tagging variable was used 
to distinguish separately the b from c and from combinatorial background events. The 



feff 

Jb s ,id 

feff 
iB s ,2D 

feff 

Jb,id 

feff 
iB,2D 



eff 
cc 
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distributions of discriminant variables are shown in Figure |8| for events selected in the sig- 
nal region, after having subtracted the corresponding distributions of background events 
obtained using events within the side-bands of the D s signal. The corresponding distri- 
butions for background events are also shown in Figure |8|. For the b-tagging variable, 
comparison with simulation is shown for the sum of the b and c events. The agreement 
between real data and simulated distributions is satisfactory. 

The use of this procedure is equivalent to increase the B° purity by 20%. 

4.4 Measurement of the Bg lifetime 

4.4.1 Proper time resolution 

For each event, the B° proper decay time was obtained from the measured decay 
length (L B o) and the estimate of the B° momentum (pb°)- The technique used is de- 
scribed in [^] . The predicted decay time distributions were obtained by convoluting the 
theoretical distributions with resolution functions evaluated from simulated events. Dif- 
ferent parameterisations were used for the two Vertex Detector configurations installed 
in 1992-1993 and 1994-1995. The proper time resolution was obtained from the (t-tj) 
distribution of the difference between the generated (£) and the reconstructed (tj) proper 
times. The following distributions were considered: 

• T^B,m(t — ti) is the resolution function for B decays with only one charmed meson 
(D s ) in the final state. It is parameterised as the sum of three Gaussian functions. 
The width of the second Gaussian is taken to be proportional to the width of the first 
one. The third Gaussian describes the component with a selected hadron coming 
from the primary vertex; the fraction (/ 3 ) of these events decreases exponentially as 
a function of the proper time: 

UbM* ~ U) = (1 - h ~ h)Gi(t - U, <n) + f 2 G 2 (t - U, a 2 ) + f 3 G 3 (t - U, a 3 ) (18) 
where G\{t — U,ai), G 2 (t — U,a 2 ) and G 3 (t — U,a 3 ) are Gaussian functions with 

o x = + ((Tpi/p) 2 x t 2 , a 2 = s X (7i, cr 3 = Jcr 2 L3 + (a p3 /p) 2 x t 2 ; f 2 and f 3 are the 

fractions of the second and third Gaussian functions, respectively; f 3 is defined as 
f 3 = exp(s 2 - s 3 t). 

The values for the decay length resolutions, aLj, the momentum resolutions, cr P j/p, 
the relative fractions, fj, and the coefficients Sj, are given in Table |6|. 

• 7^-b,2d(£ — ti) is the resolution function for B decays with two D mesons in the final 
state. In this case, the selected hadron often does not originate directly from the B 
vertex, but from the second D vertex, hence resulting in a worse resolution. This 
resolution function is parameterised in a similar way as 7^b,id(£ — t^) and the values 
of the corresponding parameters are shown in Table |6|. 

4.4.2 Likelihood fit 

The B° mean lifetime and the time distribution of the combinatorial background were 
fitted simultaneously, using selected events lying within a mass interval of ±2cr centred 
on the measured D s mass (2953 events) and events lying in the D s mass side-band (3373 
events) between 2.1 and 2.3 GeV/c 2 . The probability density function for the measured 
proper time, ti, can be written as: 

P(ti) = /^f jlD ^B s , ID (^i ) + /b12D-Pb s ,2d(^) + /b^/d-P^ID (U) + /b^D-Pb^D^) (19) 

+f: c ff p cc (u) + m g f p bkg (u). 
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D^h^ sample 


Resol. function (years) 


CTLl(ps) 




CTL3(PS) 




h 


Sl 




s 3 (ps~ x ) 


TVw(t-U) (1992-1993) 


0.149 


0.140 


0.144 


0.386 


0.15 


3.5 


-1.54 


0.14 


TVwit-ti) (1994-1995) 


0.145 


0.104 


0.169 


0.256 


0.10 


2.5 


-1.87 


0.17 


Kzoit-ti) (1992-1993) 


0.236 


0.095 


0.144 


0.386 


0.30 


3.5 


-1.21 


0.14 


TL^olt-U) (1994-1995) 


0.214 


0.094 


0.169 


0.256 


0.25 


3.5 


-1.30 


0.17 



Table 6: Fitted values of the parameters used in the resolution functions TZ B ,m(t — U) 
and TZ B ,2i)(t - U). 

The different probability densities are expressed as convolutions of the physical probabil- 
ity densities with the appropriate resolution functions: 

• For the signal: 

Pb.,w{U) = — exp(-t/r B o) <g> TZ B ,m(t - U) 

TBS 

Pb s ,w(U) = — exp(-t/r B o) <g> 1Z B ,2T>(t - U), 

where t is the true proper time. 

• For the physical background coming from all B meson decays: 

■^,id(*«) = — /B q ,iD exp(-t/r Bq )®^B,iD^-ii) 

Pb2d(U) = ^2— /B q ,2D exp(-t/r Bq ) <g> 7£b,2d(* - **)■ 

where q runs over the various b-hadron species contributing to this background and 
/B q ,iD, /B q ,2D are their corresponding fractions. 

• For the combinatorial background, the following function was used: 

PbkgiU) = f- exp(t/r-) g> G(t - U, <7_) + f+ exp(-t/r+) ® G(t - t h cr+) + 
(l-f--f+)G(t-ti,a ) 

Three distributions were considered: an exponential for poorly measured events 
having negative t (with lifetime r~), an exponential for the flying background (with 
lifetime r + ) and a central Gaussian for the non-flying one. The seven parameters 
(/~, / + , r + , r~, (T_, a + , cr ) were allowed to vary in the fit. 

• For "charm" candidates, the function P cc (ti) has the same form as Pbkg{ti) and has 
been parameterised using simulated events. In this case all the parameters were 
fixed in the fit. 

The B° lifetime fit was performed in the proper time interval between -4 ps and 12 ps 
and the result of the fit, shown in Figure |9|, is: 

r B o= 1.53 t° f 5 (stat.) ps. 

4.4.3 Systematic uncertainties on the B° lifetime 

The contributions to systematic uncertainties on the B° lifetime measurement are 
summarised in Table [7|. 
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Source of systematic uncertainty 


r B o variation (ps) 


Sample composition 


+0.013 
-0.016 


fbkg 


+0.046 
-0.050 


B s purity 


+0.005 
-0.015 


t resolution 


±0.019 


r B +(1.65±0.03 ps) 


±0.016 


r B o(1.56±0.03 ps) 


±0.014 


Analysis bias correction 


±0.040 


Total 


±0.07 



Table 7: Sources of systematic uncertainties on the B° lifetime. 

The systematic error due to the uncertainties in the relative fractions of the different 
D s sources corresponds to a ±l<x variation of the fractions / used in the likelihood fit, 
excluding f bkg which is studied separately. 

The estimate of systematics related to the evaluation of the B° purity, on an event- 
by-event basis, was obtained in the following way. The distributions of the different 
quantities contributing to the discriminant variable (Figure |8|) for signal and background 
events were re-weighted with a linear function in order to maximise the agreement between 
data and the simulation. The B s lifetime distribution was refitted with new probability 
distributions and the difference between the corresponding values of the fitted lifetime 
taken systematic error. 

Uncertainties on the determination of the resolution of the proper time receive two 
contributions: one from errors on the decay distance evaluation and the other from 
errors on the measurement of the B° momentum. The systematic error coming from 
uncertainties on the time resolution was evaluated by varying the widths ctl and o p of 
the resolution function by ±10% [|K|^7],|8|] . Finally, simulated B° events, generated with 
a lifetime of 1.6 ps and satisfying the same selection criteria as the real data, have a 
fitted lifetime of 1.64 ± 0.04 ps. The B° lifetime value obtained in the Section fl.4.2| was 
not corrected but the statistical error of this comparison (±0.04 ps) was included in the 
systematic error. The measured B° lifetime was found to be: 

r B o = 1.53±g:^(stat.) ± 0.07(syst.) ps. 

A similar analysis by the ALEPH collaboration [12|] gave a consistent result: r B o = 
1.47 ± 0.14(stat.) ± 0.08(syst.) ps. 



4.5 Lifetime difference between Bg mass eigenstates 

The B° (or B°) mesons are superpositions of the two mass eigenstates: 



|B°> = ^(|B° ) ± |B°» ; |B°) = ^( 



|B&> - |B°». 
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Neglecting CP violation, the time probability density is then given by: 

V{t) = (1 - x cp )-^- (e"™ + e-^) + x cp T L e^ (20) 

J- H + J- L 

where T L = r B o + Ar B o/2, T H = T B o - Ar B o/2. 

The first term of equation (BO) refers to final states of identified beauty flavor, as in the 



D s £ case ]I0], it does include D s +light mesons final states. The second term corresponds 
to DM+rJF- fi na l states, which are dominantly (98%) CP even eigenstates |[34|| . The 
value of x cp was taken as the ratio of -F Bs ,2D to (-F Bsi id+-^b s ,2d) : x cp = 0.22 ± 0.07 (see 
Section p|). 

Two variables are then considered: r B o (l/T B o) and Ar B o/T B o. As the statistics in 
the sample is not sufficient to fit simultaneously r B o and Ar B o/r B o, the method used 
to evaluate Ar B o /T B o consists in calculating the log-likelihood for the time distribu- 

s s 

tion measured with the Df h T sample and deriving the probability density function for 
Ar B o/r B o by constraining r B o to be equal to the lifetime (r B o = 1.56 ±0.03 ps and 



r B o/r B o = 1 ± 0(0.01) is predicted in 0). 

The log-likelihood function described in Section ^4.4.2j was modified by replacing the 
probability density for B° by (EDI) convoluted with the appropriate resolution functions. 



It was minimized in the (t b o, Ar B o/r B o) plane and the likelihood difference with respect 
to its minimum AC (Figure |l0|-a) was computed: 

ac = - \ogC tot (r B o, Ar B o/r B o) + io g £ to4 ((r B o) mi ", (Ar B o/r B ») mm ) . 

The probability density function for the variables r B o and Ar B o/r B o is proportional to: 

P(r B o,Ar B o/r B o) cxe" A£ . 

The Ar B o /r B o probability distribution was then obtained by convoluting P(r B o , Ar B o /r B o) 
with the probability density function /( T 0=T )(t b o), expressing the constraint r B o = r B o, 

B s B d d 

and normalising the result to unity: 

/ V (r B o , Ar B o /r B o ) / (T =T ) (r B o ) dr B o 

p(Ar B o/r B o) 



/ V (r B o , Ar B o /r B o ) f (T =T ) (r B o ) dr B o dAT B o /T B o 

■ d 

where 

/(r B o =r B o ) ( r B0 ) = 1 / ( V^^o ) exp ( - (r B o - T B o ) 2 / 2(7 2 ) . 

s d d d B d 

The upper limit on Ar B o/r B o, calculated from P(Ar B o/r B o), is: 

Ar B o/r B o < 0.67 at the 95% C.L. 

This limit takes into account both statistical uncertainties and the systematic coming 
from the uncertainty^ on the B° lifetime. 

The systematic uncertainty originating from other sources was evaluated by convo- 
luting the probability function V(t b o, Ar B o/r B o) with the probability function of the 
corresponding parameters: V(AT B o/T B o) = 

J V(r B o , AT B o /r B o , xl ys , . . . , x n sys ) f\ TB0 =Tb0 ) (r B o ) / {xl ys )... f\x™ ys ) g?t b o dx\ ya . . . dx n gys 

s d 

/ V (r B o , Ar B o /r B o , x\ ys , . . . , x n sy s ) / (T =Tb0 ) (r B o ) / {x x sy s ) . . . f { x n sy s ) dT B odx 1 sys ...dx n sys dAT BO jT B o 

s d 

5 The uncertainty due to the theoretical prediction of the equality of the r B o and r B o lifetimes is negligible with respect 
to the present error on r B o . 

d 



23 



where x l sys are the n parameters considered in the systematic uncertainty and f{x l sys ) are 
the corresponding probability densities. 

Only three systematics were considered here: the relative fraction of combinatorial 
background fbkg, the B° purity of the selected sample and the x cp fraction. Other sys- 
tematic uncertainties are expected to be small as they are in the lifetime measurement. 

The Ar B o /r B o probability distribution, obtained with the inclusion of the systematics, 
is shown in Figure |TU]-c, the most probable value for Ar B o /r B o is 0.35 and the upper limit 
at the 95% confidence level is: 

Ar B o/r B o < 0.69 at the 95% C.L. 

It should be noted that the world average of the B° lifetime cannot be used as a constraint 
in such an analysis, since it depends on Ar B o and on r B o. Moreover, this dependence is 
also different for different decay channels. In the D^h T case the expression of the average 
B° lifetime is given by: 

r B .(Df h*) = (1 - *„) \^^ff m + Wr L • (21) 



4.6 Study of B°-B° oscillations 
4.6.1 Tagging procedure 



The tagging algorithm was explained in Section |3.5.1| . An event is classified as a mixed 
or an unmixed candidate according to the relative signs of the D s electric charge, Qd, 
and of the tagging purity variable, x tag - 

Mixed candidates have x tag -Qd < 0, and unmixed ones x tag • Qd > 0. The probability, 
ej^ 9 , of tagging the b or the b quark correctly from the measurement of Xt ag was evaluated 
using a dedicated simulated event sample. The average tagging purity of the x tag variable, 
given by the simulation for true Bj? — > D~h + X decays, is (71.4±0.4)%. The purity is lower 



than that obtained in the D s l sample fl(| because not all B s charged decay products are 
reconstructed in the present analysis. It was verified that the tagging purity is the same 
for different B hadron species and varies by less than about ±2% whether the B° has 
oscillated or not. This effect is taken into account in the systematics. The corresponding 
probability distribution for events in the combinatorial background was obtained using 
data candidates selected in the side-bands of the D s signal: the probabilities of classifying 
these events as mixed or as unmixed candidates are called e^g and e^p tx , respectively. 
For the charm events, the analogous probabilities are called e™ lx , e^ mix and their values 
were taken from the simulation. 



4.6.2 Fitting procedure 

From the expected proper time distributions and the tagging probabilities, the prob- 
ability functions for mixed and unmixed event candidates are 

pmix(± \ _ feff pmix (± \ , feff pmix (+ \ j_ feff r>mix (+ \ j_ feff prnix l + \ (00\ 
r \ L i) — JB B) m r B s ,lD\ L i) ~r iB s ,2D- r B s ,2Dl''«J ~>~ J B,lT> r B,lD\ L i) ~>~ JB,2D- r B,2DV t V 

_l_ feff pmix ( j. \ _i_ feff pmix ( + \ 
"•"ice -'cc "r Jbk g r bk g \ L %)i 

where ti is the reconstructed proper time. The analytical expressions for the different 
probability densities are given in the following, with t being the true proper time: 
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• B° signal mixing probability: 

P£%(ti) ={ ef 9 V^(t) + (l-ef 9 )V^(t) }®K B , m {t-U) (23) 

• Physical background mixing probabilities: 

*b&(*<) = {/B Sl 2D4 a 7^o exp(-t/r B o)} ® ft B , 2D (t - *<) (24) 

^bTd(^) = { /B dl iD( 6^P^ 1D (t) + (1 - ef 9 )7~(0 ) + 

/b+,id(1 - ep/r B+ exp(-t/r B +) + (25) 
/A b ,iD(l - 4° S )/ r A b exp(-t/r A J } <g> 7^b,id(* - *i) 

^bTd(^) = { /B d)2D ( ei a9 Ppi?(t) + (1 - ef s )P™* 2D (t) ) + 

/b+,2d(1 - 4 ag )/ T B+ exp(-t/r B +)+ (26) 
/A b ,2D(l - 4 a9 )/ r A b exp(-t/r A J } <g> 7£ B)2 d(£ - *») 

• Mixing probability for charm component: 

= c^cc(^) (27) 

• Combinatorial background mixing probability: 

PbkTiU) = £bkg'Pbkg(ti) (28) 

The oscillation analysis was performed in the framework of the amplitude method |13| 
as described in Section |3.5.3| . Considering only statistical uncertainties, the limit is: 

Am B o > 4.2 ps" 1 at the 95% C.L. (29) 

with a corresponding sensitivity equal to 3.1 ps _1 . At Am B o = 10 ps _1 , the error on the 



amplitude is 2.3 (see Figure |Tl"l) . 



4.6.3 Study of systematic uncertainties 

Systematic uncertainties were evaluated by varying the parameters which were kept 
constant in the fit, according to their measured or expected errors. 

• Systematics from the tagging probability: 

a conservative variation of ±3% on the expected tagging probability for the signal 
and for the other three processes contributing to the D^h^ candidates was used. The 
same variation is assumed for the tagging purity for the charm and combinatorial 
background events. The central values of these purities were fixed to the simulated 
ones. 

• Systematics from the B° purity: 

the same procedure already applied to the lifetime measurement has been used. 

• Systematics from the resolution on the B decay proper time: 

the same procedure already applied to the lifetime measurement was used. In ad- 
dition, the systematic error due to the variation of the proper time distribution of 
the combinatorial background was considered: the parameters used to define the 
background shape in the lifetime fit were varied according to their fitted errors. 

The inclusion of systematic uncertainties lowers the sensitivity to 2.7 ps" 1 and the 
95% C.L limit becomes Am B o > 4.1 ps _1 . 



The analogous analysis has been performed by the ALEPH collaboration |T2[], which 



set a limit at Am B o > 3.9 ps 1 at the 95% C.L., with a better sensitivity of 4.1 ps 1 . 
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5 Combined limit on Am B o using exclusive Bg and 
D^h T events 

This paper presents two analyses on Am-so using exclusively reconstructed B° mesons 
and D^h^ events. Their results have been combined (Figure |12|), taking into account 
correlations between systematic uncertainties in the two amplitude measurements J7J. A 
limit at the 95% confidence level is obtained: 

Am B o > 4.0 ps" 1 at the 95% C.L. (30) 

with a corresponding sensitivity equal to 3.2 ps -1 (with statistical errors only, the limit 
would be AmBO > 4.0 ps -1 at the 95% C.L. with a sensitivity of 4.4 ps -1 ). Figure |12 
also shows the error on the amplitude for different values of Am B o. 



6 Combination of the DELPHI Am B o, r B o and 
Ar B o/r B o analyses 

DELPHI has performed three analyses on Am B o using candidates [JjlJ, D^h T 

events and exclusively reconstructed B° mesons. They were combined, taking into ac- 
count correlations between the event samples and between systematic uncertainties in the 
different amplitude measurements (Figure |i~3|a). This gives the following limit for Am-go: 

Am B o > 4.9 ps -1 at the 95% C.L. 
with a sensitivity of Amso = 8.7 ps -1 

The exclusion probability for this limit is 88%. The variation, with Ameo of the error on 
the amplitude is given in Figure |T3]b. 

The results of two analyses on r B o and Ar B o /r B o using Df£^ and D^h T events provide 
the following DELPHI results: 

r B o = 1.46 ±0.11 ps 
Ar B o/r B o < 0.45 at the 95% C.L. 



7 Conclusion 

Using about 3.5 million hadronic Z decays registered by DELPHI between 1992 and 
1995, two samples of events have been selected. The first one consists of 44 reconstructed 
B° events: 11 candidates (including 30% of background) are completely reconstructed 
and 33 candidates (including 55% of background) are partially reconstructed (7r° and/or 
7 are not detected). This analysis used twelve different decay channels of the B° meson 
and is a first attempt to use such events for the oscillation studies. Due to the excellent 
proper time resolution, this sample gives some contribution in the high Am^o region. 

The second sample contains 2953 D^h T candidates (including 60% of background) 
with completely reconstructed D~ mesons in the 07r~ and K*°K~ decay channels. Using 
the D^h T sample, three studies have been performed. The B° lifetime has been measured 
and a limit on the fractional width difference between the two physical B° states has been 
obtained: 

r B o = 1.53l£^(stat.) ± 0.07(syst.) ps 
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Ar B o/r B o < 0.69 at the 95% C.L. 

This last result has been obtained under the hypothesis that r B o = r B o. 

Combining the two studies on B° — B° oscillations, a limit at the 95% C.L. on the 
mass difference between the physical B° states has been set: 

Am B o > 4.0 ps" 1 at the 95% C.L. (31) 

with a corresponding sensitivity equal to 3.2 ps _1 . 

Previous DELPHI results on B° lifetime obtained with the D^h^ sample 
superseded by the analysis presented in this paper. 

Combination of the DELPHI Am B o , r B o and Ar B o /T B o analyses gives: 

Am B o > 4.9 ps" 1 at the 95% C.L. 
with a sensitivity of Am B o = 8.7 ps _1 

r B o = 1.46 ±0.11 ps 
Ar B o/r B o < 0.45 at the 95% C.L. 
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Figure 1: B° mass spectrum for the candidates selected in the twelve decay channels 
described in Section \3.0\ . The data are indicated by the points with error bars and the 
result of the fit has been superimposed. Details on this fit are given in the text. The 
histograms represent the expected contribution from beauty events (after having removed 
the exclusively reconstructed B° decay channels), from charm events and from light quark 
events. 
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Figure 2: The Monte Carlo composition of the satellite peak including the effects of 
experimental resolution. The notations in the plot are the following: D s p corresponds to 
Dj7p + and D* _ p + decay channels of the B°; D s 7r(a!) corresponds to D*~7r + and D*~a^/ 
D°K7r(ai) shows the contribution from D (2007)°K~7r + and D (2007)°K~a^ decay chan- 
nels. All contributions have been normalised according to the evaluation of the branching 
fractions discussed in Section and including reconstruction efficiencies evaluated using 
simulated events. 
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Figure 3: Mass spectra for a) B+ -> D°7r+ ; b) B° d -> D*(2010)-tt+ and 
Bj} -> D*(2010)-a^ decays. In the first plot a signal from B+ -> D*(2007)°7r+ ; 

B + — > D p + and B + — > D*(2007)°p + (satellite peak) is also visible. Details on the fit 
are given in the text. 
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Figure 4: B° mass spectrum for the sum of B[J — > D n~n + and B° — > D 7r~ai decays 
selected in the four decay channels described in Section TTie data are indicated by 
the points with error bars and the result of the fit has been superimposed. The histograms 
represent the expected contribution from beauty events (after having removed the exclu- 
sively reconstructed decay channels), from charm events and from light quark events. 
The widths of the signals have been fixed according to the values found in the simulation. 
Details on the fit are given in the text. 
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Figure 5: Exclusive B° analysis: variation of the oscillation amplitude A as a function 
of Am^o. The lower continuous line corresponds to A + 1.645 o~(A) where cr(A) includes 
statistical uncertainties only, while the shaded area shows the contribution from system- 
atics. The dashed-dotted line corresponds to the sensitivity curve. The lines at A=0 and 
A=l are also given. The points with error bars are real data. 




Figure 6: °~{A) as a function of Am B o. a) Exclusive B° analysis: the full curve 
shows the result from the data. The dashed curve shows the average result from 100 toy 
experiments with the same statistics as in data and the shaded area gives the ±2a(A) 
region around this average. The systematic effects are not included, b) Comparison of 
the o~(A) as function of Atti-qo for three analyses: the full curve shows the result from the 
exclusive B° analysis, the dashed and dotted curves show the result from D^h T and D^£ T 
analyses, respectively. 
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Figure 7: Invariant mass distributions for D s candidates selected in (\>ix~ and K*°K~ 
decay channels. The upper and lower plots refer to data samples registered in 1994-1995 
and 1992-1993, respectively. The selected D s candidates are accompanied by a hadron of 
opposite charge (or by several hadrons), measured in the same event hemisphere. The 
curves show the fits described in the text. 
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Figure 8: Distributions of the variables used to increase the B s purity. The distribu- 
tions on the left are for events selected in the signal region after having subtracted the 
corresponding distributions of background events, which have been obtained using events 
situated in the side-bands of the D s signal. The corresponding distributions for background 
events are shown on the right. The points with error bars correspond to the data and the 
histograms are simulated events. For the | cos -0 1 distribution, only the 07r + decay mode 
is shown because the cut on this variable was set at 0.6 for K*°K~ and at 0.4 for <pii + 
channel. In addition, the D s mass was used as the fifth discriminating variable: the signal 
has a Gaussian and the background has an exponential distribution. 
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Figure 9: D^h T analysis. Upper plot: Proper time distribution for events in the signal 
mass region. The points show the data and the shaded regions correspond to the different 
contributions to the selected events. The curve shows the result of the fit described in 
the text. Lower plot: the same as the upper plot but for events situated in the D s mass 
side-band. 
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Figure 10: h T analysis: 68%, 95%, 99% C.L. contours of the negative log-likelihood 
in the plane Ar B o/r B o — r B o a) without and b) with r B o = r B o constraint. The point 
indicates the minimum, c) Probability density distribution for Ar Bs /r Bs ; the two lightly 
shaded regions at the 68% C.L. and dark one at the 95% C.L. are also shown. 
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Figure 11: D^h T analysis. Left plot: variation of the oscillation amplitude A as a function 
o/AmBO. The lower continuous line corresponds to A + 1.645 o~(A) where o~(A) includes 
statistical uncertainties only, while the shaded area shows the contribution from system- 
atics. The dashed-dotted line corresponds to the sensitivity curve. The lines at A=0 and 
A=l are also given. The points with error bars are real data. Right plot: variation of the 
error on the amplitude as a function of Ameo , including systematic uncertainties. 
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Figure 12: Combination of the D^h T and exclusive B° analyses. Left plot: variation of the 
oscillation amplitude A as a function of Am B o . The lower continuous line corresponds 
to A + 1.645 cr(A) where o~(A) includes statistical uncertainties only, while the shaded 
area shows the contribution from systematics. The dashed-dotted line corresponds to the 
sensitivity curve. The lines at A=0 and A=l are also given. The points with error bars 
are real data. Right plot: variation of the error on the amplitude as a function of Am B o 
including systematic uncertainties. It should be noted that the D^h^ analysis dominates 
at low values of Am B o and the exclusive B° analysis dominates at large values of Am B o . 
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Figure 13: Combined DELPHI analysis: a) variation of the oscillation amplitude A as 
a function of Atobo. The lower continuous line corresponds to A + 1.645 o~(A) where 
o~(A) includes statistical uncertainties only, while the shaded area shows the contribution 
from systematics. The dashed-dotted line corresponds to the sensitivity curve. The lines 
at A=0 and A=l are also given. The points with error bars are real data, b) variation 
of the error on the amplitude as a function of Atobo . 



